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Abstract 

We show that the thermodynamic Bethe ansatz equations for one-dimensional 
integrable many-body systems can be reinterpreted in such a way that they only 
code the statistical interactions, in the sense of Haldane, between particles of 
identical or different momenta. Thus, the thermodynamic properties of these 
systems can be characterized by the generalized ideal gases recently proposed 
by one of us. For example, the Yang- Yang <5-function gas is a gas with spe- 
cific statistical interactions between particles of different momenta, while the 
Calogero-Sutherland system provides a model for an ideal gas of particles with 
a fractional statistics. 
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It is well-known that in 1 space dimension (or in 1 + 1 space-time dimensions), 
bosonic theories can be mapped into fermionic ones and reciprocally [I], |2|]. Further- 
more, some Id many-body systems even exhibit a continuous boson-fermion inter- 
polation, when the coupling constant varies in appropriate range. This is known to 
occur when the system can be solved by the thermodynamic Bethe ansatz, such as 
the Yang- Yang S- function gas [|[] or the Calogero-Sutherland system ||J]. In this note, 
we give a simple characterization of the thermodynamic properties of these systems 
in terms of Haldane's statistical interactions || and the consequential generalized 
ideal gases ||. In particular, we point out that the Bethe ansatz equations in these 
models can be re-written in such a way that all the dynamical interactions have been 
transmuted into statistical interactions. 

1 Fractional statistics and thermodynamics 

Haldane || has recently introduced a notion of fractional statistics which is inde- 
pendent of the dimension of space. This notion is not based on the monodromy 
properties of the iV-particle wave functions, but on the way the number of available 
single-particle states varies when particles are added into the system. More precisely, 
consider a system with a total number of particles N = J2j Nj, with Nj the number 
of particles of the species j. Consider now adding a particle of the species i into 
the system without changing its size and the boundary conditions. Keeping fixed 
the positions of the N particles of the original system, the wave function of the new 
(N + l)-body system can be expanded in a basis of wave functions for the added 
particle. We denote by Di the dimension of this basis. The important point is that 
this dimension may depend on the numbers Nj of particles in the original system. 
Assuming that this dependence is linear, Haldane defines the 11 statistical interaction! 1 
through the relation || : 



Clearly, for bosons the numbers of available single-particle states are independent of 
the numbers Nj, and g^j = 0. For fermions, the numbers of available single-particle 
states decrease by one for each particle added, and = 6y. 



dDi 




dNj 



~9ij ■ 
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One of the ideas underlying the introduction of the generalized Pauli principle (|I|) 
is the fact that bosons and fermions can be considered on an equal footing as far as 
state counting is concerned. Indeed, for bosons or fermions, the number of states of 
(N + 1) identical particles distributed among G accessible orbitals is : 

W - ^/ + iV ) ! 

bJ (N + l)\(D bJ -l)\ 

with 

D b (N) = G, and D f (N) = G - N . 

One of us || has recently generalized this to fractional statistics by directly count- 
ing the many-body states and assuming that the total number of states with {Nj} 
particles is : 

V m [D.mri - 1]! • 1 > 

where Di({Nj}) is obtained by integrating ([]]) : 

A({iv,}) + E^^ = ^°, (3) 

3 

with G° = Di({0}) a constant, which is interpreted as the number of available single- 
particle states when no particle is present in the system. Namely, G° are the bare 
numbers of single-particle states. In ref. 0, for % ^ j was termed as "mutual 
statistics" . 

Knowing how to enumerate states, it is then possible to study the thermody- 
namics. In the thermodynamic limit, the numbers of particles {Nj}, as well as the 
bare numbers of available states {G®}, become infinite. But the occupation numbers 
rij = (Ni/G®) remain finite. The entropy is S = kB^ogW with ks the Boltzman 
constant. By definition, we call a system a generalized ideal gas, if its total energy 
with (Nj) particles is simply given by 

£ = £M (4) 
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with constant e°. For such gases, the thermodynamic potential f2 = —PV at equilib- 
rium can be evaluated by minimizing 

n = E - TS - J2 N jM 

j 

with respect to the variation of the densities rij. Here T is the temperature and //, 
the chemical potential for the species i. 

The resulting thermodynamics can be summarized as follows |J : 

Q = -k B Tj2G° l lo g ( 1 -±^-) , (5) 

i \ Wi J 

where the functions w% are determined by the equations : 

log (1 + v,i) + £ 9H ^ [j^-] = 4^ • ( 6 ) 
j \ 1 + Wj J k B T 

These relations completely specify the thermodynamics of the generalized ideal gas. 
The other thermodynamic quantities can be derived from the relation : 

dVl = -SdT - £ Nidfn - PdV . 

i 

In particular, the occupation numbers n« are obtained from riiG® = —dVt/d^i. It 
gives : 

"■EM, • < 7 > 

3 

where B is a matrix with entries : Bij = Wi5ij + hij, with G®hij = gijG®. 

The simplest example considered in is the ideal gas of particles with a fractional 
statistics; i.e. g^ = gS^j, = /i. In this case, the statistical distribution of rij is then 
given by n, = l/(wi + g), with Wi satisfying eq.(|6]) which now becomes: 

Wi 9 (1 + w^ 1 - 9 = exp ((e? - fi)/k B T) . (8) 

For g — (or g — 1), we recover the bosonic (or fermionic) occupation numbers. 
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2 The Bethe ansatz equations and state counting 



The thermodynamic Bethe ansatz (TBA) is a method developed by Yang and Yang 
H for finding the thermodynamic properties of a one-dimensional (Id) gas of particles 
whose interaction is modeled by an integrable hamiltonian. For simplicity, we will 
assume that the particles have no internal quantum number. 

In the Bethe ansatz approach to integrable models, the information is encoded in 
the two-body S-matrix. We denote it by S(k) with k the relative momentum of the 
scattering particles; we have S(k) = — exp(—i6(k)), where 8(k) is the phase shift and 
it is odd in k. The eigenstates of the iV-body hamiltonian (with the periodic boundary 
condition) are labeled by N (pseudo-)momenta {k r } (r = 1, , ■ • • , N) (defined in the 
asymptotic regions), which are solutions of the Bethe ansatz equations : 

e ik r L = -Q S ( ks _ ^ for aU r ^ (g) 

with L the length of the system. The energy of the eigenstate \k±, ■ • • , k^) is found 
to be : 

£ = E*, (io) 

r 

where e°(k) is some universal function, e.g. e°(k) = k 2 . 

Using this ansatz, Yang and Yang have derived the thermodynamic potential of 
such gas, which we denote by Q = Q(L, T, /i), at temperature T and chemical potential 

j = -k B T[°° ^log(l + e -'W/*^ , (11) 
where the function e(k) is determined by solving the TBA equations : 

e(*) - k B T r ^ 4>{k, k') log (l + e-< k '^ T ) = e°(k) - // , (12) 
with <j)(k, k') the derivative of the phase shift : 

0(fc, k') = i— log S(k, k') = 6'(k - k') . (13) 
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Here we point out that a simple characterization of the statistical properties of 
these systems can be given in terms of the notions we reviewed in the last section. 
Namely, the thermodynamic properties of the Id TBA gases can be characterized as 
generalized ideal gases with non-trivial statistical interactions (or mutual statistics), 
if the particles with different momenta are considered as belonging to different species 
(subject to one and the same chemical potential). 

As we now explain, the information about statistical interactions is encoded just 
in the Bethe ansatz equations (^), if it is rewritten in approriate form. Indeed, taking 
as usual the logarithm of eq.(§), the Bethe ansatz equations become : 

Lk r = ^2 6{k r — k s ) + 2n I r , (14) 

with 8(k) the phase shift. Here {/ r } (r = 1, • • • , N) is a set of integers or half integers, 
depending on N being odd or even, which one may choose to serve as the quantum 
numbers labeling the eigenstates, instead of the momenta {k r }. For simplicity, let us 
consider the cases when these 7 r 's can be chosen to be all different and still provide a 
complete set of solutions (i.e., one has a fermion description for the purpose of state 
counting). Though in this description the statistics looks simple (just like fermions), 
the dynamics is nontrivial in that the total energy is a complicated function of {I r }- 



Now let rewrite eq. (|T4 ) , in the thermodynamic limit (with iV — > oo, L — ► oo and 
N/L fixed), in the momentum description for state counting, by dividing the range 
of momentum k into intervals with equal size Ak and labeling each interval by its 
midpoint fej. (At the end of the discussion, we will take the limit Ak — > 0.) Then 
we treat the particles with momentum in the i-th interval as belonging to the i-th 
species, for which the bare available single-particle states is obviously G° = LAk/2it. 
Now let us introduce the distribution for the roots, k r , of the Bethe ansatz equations 
©by: 

Lp(ki)Ak = Ni = No. of k r in - Ak/2, k { + Ak/2) . (15) 

Note that 2np(ki) = Ni/G® is the occupation number distribution. 
Now, motivated by eq. (0), let us define : 

— J(Jfe) =k--Y,6(k-k B ) = k-Y l 6(k- k^pik^Ak . (16) 

k 3 j 
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It is easy to see that whenever I{k) = I with / an integer, the corresponding value of k 



is a possible root of the Bethe ansatz equations fll4|) , and thus represents an accessible 
orbital. Therefore, the number of accessible orbitals, Di({Nj}) = L p t (ki)Ak, in the 
i-th momentum interval is given by 

Di = I(ki + Ak/2) - I(ki - Ak/2) ; (17) 

or, using eq. ([16]) and taking the Ak — > 0, we have : 

1 1 r°° 

Pt{k) = 2^ " 27 loo dk 6 ' {k ~~ k ' )p{k,) ■ (18) 
Recall that in the fermion description for state counting, the same W in eq.(f2|) is 
obtained with the number of accessible orbitals for the i-th species taken to be 

Di({Nj}) = DidNj}) +N i -l = G° i +N i -l-J29i j N 3 . (19) 

j 

Comparing eq. (jTST) with the continuum form of this equation, we identify the statis- 
tical interactions in the momentum description to be 

g(k,k') = 6(k-k') + —6'(k-k'). (20) 

This shows that the dynamical interaction, which is summarized in the two-body 
phase shift, is transmuted into a statistical interaction in the momentum description. 

The statistical interaction between particles of identical and different momenta is 
a manifestation of the following features of the spectrum of a TBA gas: The ground 
state corresponds to an equidistribution of the integral quantum numbers {/ r } in a 
certain interval: I r+ i — I r — 1; and the excited states correspond to particle/hole 
excitations in the integral lattice for /: 

I r+1 -/ r = l + M r \ (21) 

where is the number of holes, i.e. unoccupied integer numbers, between J r+1 and 
I r . Therefore, changing to the momentum description, we have the total density of 
accessible states 



p t (k) = p(k)+p h (k) , 
7 



(22) 



where p(k) is the particle density defined in eq.(|15D, and Ph(k) = M h {k)p{k) is the 
hole density, which is related to the distribution of unoccupied sites in the integer 
lattice of quantum number /. This equation is essentially the Bethe ansatz equation 
in the thermodynamic limit. In the light of this, eq.(|l8|) implies the linear dependence 
of the hole density Ph{k) on the particle density p(k'). 

Having identified gij or g(k, k'), we turn our attention to the condition ([|) for a 
generalized ideal gas. Because of eq. (PH|), this condition holds good: 

- = / dk p{k)t°{k) . 

LJ J — OO 

Thus there is no interaction energy between particles of different momenta. Contrary 
to the fermion description, in the momentum description the dynamics looks simple 
but statistics is non-trivial. 

Now we can apply the formulas for the generalized ideal gas in the last section to 
derive the thermodynamic properties of the TBA system. It is not surprising that 
they lead to exactly the same results as Yang- Yang's. For example, eqs. @ and (||) 
coincide with eqs. (|Tl~r) and (13), with the following idenfication: 



GVL — p°{k) = -|- (23) 

Z7T 

Wl — Y(A0 = exp(j^) (24) 

g .. <__> g{k,hf) = 5{k-hf) + ^{k,k') (25) 

with the discrete sum replaced by the integral over momentum. In particular, Yang- 
Yang's form of the thermodynamic limit of the Bethe ansatz equations : 

'•+ 00 dk! 



Ph 



^<P(k,k')p(k')=p°(k) = - (26) 

-oo Z7T Z7T 



agrees with our state-counting equation ([19]) or (Tl8|), while their entropy : 
S f+°° 

-= dk [(p + p h ) log(p + p h ) -p log p-p h log p h ] (27) 

Li J —oo 

is just 5 = kB^ogW, with W given by our eq.(0), with the correspondence : 

Nt/L <— > p(k) (28) 
Di({Nj})/L > Ph (k) . (29) 



This last idenfication is quite natural: the hole density ph(k) clearly represents the 
density of available states for an additional particle to be added. 



3 Examples of generalized ideal gas 

The simplest example of the TBA system is the repulsive ^-function bosonQ gas, with 
a two-body potential: V(x) = 2c5(x), c > (x being the relative coordinate). This 



model is integrable with a fermion description for the Bethe ansatz equations (14) 
with the two-body phase shift : 

6{k) = -2tan'\k/c) . (30) 

Therefore, <fi(k, k f ) = —2c/[(k — k!) 2 + c 2 ], and the statistical interaction is : 

»(*.*0 = f(*-v)- ^ + (t e _ y) , . (3D 

If c — > 0, the second term becomes —S(k — k'), so we have an ideal boson gas, as it 
should be. If c — > oo, the second term vanishes for all finite k, and we have an ideal 
Fermi gas: g(k,k') = 8(k — k'). For intermediate values of c, the thermodynamics 
of the gas is that of a generalized ideal gas with specific mutual statistics between 
different momenta. 

Another well-known example is provided by the Calogero- Sutherland model 
This is a model of particles interacting through a 1/r 2 potential. It is integrable, and 
its two-body S-matrix is 

S(k) = - exp [-z7r(A - 1) sgn{k)\ , (32) 

with A the coupling constant. Therefore, <fi(k,k') = 2it(\ — l)8(k — k'), and the 
statistical interaction is : 

g(k, k') = X5(k - k') . (33) 
The Bethe ansatz equation then reads : 

= p h (k) + Xp(k) . (34) 

Z7l 



^^Boson here specifies only the symmetry of many-body wave functions, and does not refer to the 
statistical interaction that determines or characterizes the thermodynamic properties of the gas. 
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The coupling constant A also governs the statistical interaction : if the density of 
particles of momentum k increases by a unit, then the holes density decreases by A. 
The bare energy is e°(k) = k 2 , and the TBA equations reduce to algebraic equations 

Y(k) x (1 + Y(k)) 1 ~ x = exp ((e°(A;) - n)/k B T) . 

Eq.fl33|) shows that the statistical interaction is purely between particles with 
identical momenta. In this respect, the Calogero-Sutherland system appears clearly 
as an ideal gas of particles with a fractional statistics^, exhibiting the phenomenon 
of fractional exclusion || . 

Besides the values A = 0, which corresponds to free bosons, the values A = 
1/2, 1, 2 are special, since at these values this system is related to the system 
of orthogonal, unitary, symplectic random matrices respectively. The case A = 1 
corresponds to free fermions. The case A = 1/2 corresponds to "semions" with a 
half statistics. The case A = 2 is dual to the semionic case: 
Y^ x=2 \k\ — T) e( e °~^/ fesT . In these cases, the TBA equations can be simply solved. 
Finally, from eq.(|3^) we see that the duality A < — > 1/A ||, which exchanges the 
system with coupling A and 1/A, corresponds to exchanging particles and holes. 

4 Possible generalizations 

First let us remark that in some models a boson description for the Bethe ansatz 
equations is possible: i.e. the integers I r in eq . ([14]) are allowed to repeat. Our 
discussions can be applied to these cases. Then the first term in eq. (|^) disappears. 

There are at least two potential generalizations of this approach: 
- We can generalize this approach by relaxing the assumption we made in eq. (Q), 
i.e. the energy spectrum for each particle species is degenerate. Let us denote by 
e°(aj) the energy spectrum with ctj = 1, • • ■ then an alternative anstaz for the 
2 The relation between fractional statistics and this model has been discussed, e.g., in ref. 0. 



10 



energy would be : 

J5=E^i(«i)^(«i) . 

where Nj(atj) are the occupation numbers with J2 aj ^j( a j) = iVj-This generalization 
should be related to the non-scalar Bethe ansatz (for models with internal degrees of 
freedom) . 

- We can use this approach to generalize the Fermi liquid theory by incorporating 
statistical as well as dynamical interactions. At thermodynamic equilibirum, the 
system will be described by the particle density p(k) and the hole density ph{k). In 
D space dimension, the entropy will be assumed to be the same as in eq.(P?D : 

S 
L 

where the densities p(k) and ph{k) are assumed to be coupled with a phenomenological 
statistical interaction g(k, q): 



J d D k [{p + p h ) log(p + p h ) - p log p - p h log p h ] , (35) 



p h (k)+ d D qg(k,q)p(q)=p° . (36) 



Mimicing the Fermi liquid theory, the simplest ansatz for the energy is quadratic in 

the particle density : 

dE 
dp(k) 

where V(k,q) is a phenomenological dynamical interaction. As we will report else- 
where, the thermodynamic properties can be computed from eqs. (|35|, |36| , |37j ). How- 



e°{k) + J d D qV(k,q)p(q) = e dr (k; p) , (37) 



ever, since the fields creating the particles with fractional statistics may possess frac- 
tional anomalous dimensions, we cannot guarantee that the ansatz ( |3"7D is adapted to 
the description of the low-energy/long-distance behavior of a system of such particles. 
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Note added: After the completion of the paper, we learn that since last year, the 
problem of state counting from BAE for the non-scalar integrable models (for particles 
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with internal degrees of freedom) has been studied in several papers, in some of which 
the connection to the fractional statistics in Haldane's sense was noted, though no 
thermodynamics was addressed. See ref. 0. We thank B. McCoy for telling us the 
existence of these papers. 

References 

[1] D.C. Mattis and E.H. Lieb, J. Math. Phys. 6 (1965) 304. 

[2] S. Coleman, Phys. Rev. Dll (1975) 2088. 

[3] C.N. Yang and CP Yang, J. Math. Phys. 10 (969) 1115. 

[4] F. Calogero, J. Math. Phys. 10 (1969) 2191, 2197; 
B. Sutherland, J. Math. Phys. 12 (1971) 251. 

[5] F.D.M. Haldane, Phys. Rev. Lett. 67 (1991) 937. 

[6] Y.S. Wu, University of Utah preprint, February 1994. 

[7] See, e.g., F.D.M. Haldane and B.S. Shastry, in the proceedings of the 16th 
Taniguchi Symposium, eds. N. Kawakimi and A. Okiji, Springer Verlag 1994, and 
references therein. 

See also, P. Polychronakos, Nucl. Phys. B234 (1989) 597. 

[8] R.P.Stanley, Adv. Math. 77 (1989) 76; 
M. Gaudin, Sacaly preprint SPhT-92-158. 

[9] See, for example, S. Dasmahapatra, R. Kedem, B. McCoy and E. Melzer, J. Stat. 
Phys. 74 (1994) 239; 

R. Kedem and B. McCoy, Stony Brook preprint; 

S. Dasmahapatra, R. Kedem, T.R. Klassen, B.M. McCoy, and E. Melzer Int. J. 

Mod. Phys. B 7 (1993) 3617; 

E. Melzer, preprint hepth 9312043; 

A. Berkovich, preprint hepth 9403073. 



12 



